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We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of 
the color group, and present the implementation of the RHMC algorithm for simulating dynamical 
Wilson fermions. A first dataset is presented for the SU(2) gauge theory with two fermions in the 
adjoint representation, which has been proposed as a possible technicolor candidate. Simulations 
are performed on 8^ x 16 lattices, at fixed lattice spacing. The PCAC mass, the pseudoscalar, 
vector and axial meson masses, the pseudoscalar meson decay constant are computed. The ex- 
trapolation to the chiral limit is discussed. However more extensive investigations are needed in 
order to control the systematic errors in the numerical results, and then understand in detail the 
phase structure of these theories. 
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1. Introduction 

Gauge theories with fermions in two-index representations play important roles in different 
scenarios for Physics beyond the Standard Model (SM): supersymmetric theories, technicolor mod- 
els [|l|, ^,|3|,|^, orientifold theories [|5|,^. So far the lattice is the unique tool for studying strongly- 
interacting theories from first principles. Moreover recent algorithmic progress and the increase in 
computer performance allow these theories to be studied by Monte Carlo simulations at an afford- 
able cost. On the other hand the expertise developed in many years of simulating lattice QCD must 
be exported with some care as we try to interpret results for generic strongly-interacting theories, 
since these may reveal a very different phenomenology from QCD. 

Therefore the challenge this work attempts to address is twofold. We generalize the standard 
RHMC algorithm for simulating QCD to the case of SU(A'^) gauge group and fermions in a generic 
representation. We briefly describe the main characteristics of our bespoke code - HiRep - in 
which we have implemented the two-index representations for generic number of colors. Then we 
present some preliminary result for SU(2) with two Dirac fermions in the adjoint representation. 

The number of fermions in this theory is close to the estimated critical value where an infrared 
(IR) fixed point is expected to appear [^. If an IR fixed point exists, the theory is conformal in the 
IR, in particular it does not confine and does not break chiral symmetry. If the theory does not have 
the IR fixed point, the vicinity to the critical number of flavours can induce a walking behaviour 
of the renormalized coupling constant, which would make such a theory a possible realization 
of the technicolor mechanism [||, ^]. Which scenario is the correct one can only be established 
by numerical simulations on the lattice. However for realistic numerical simulations, the finite 
volume, the non-zero fermion mass, and the lattice artifacts introduce systematic errors that are 
likely to obscure the behaviour one would naively infer from the existence of an IR fixed point. In 
order to obtain clear and robust results, we need to go to large volumes, small masses and small 
lattice spacing, keeping in mind that we do not know a priori what large and small mean in this 
case. 

Presenting the results, we will try to emphasize which are compatible with a standard QCD- 
like scenario, and which deviate from it giving hints of new interesting dynamics. 



At the end of this introduction, we want to point out that the same theory was explored in [ |10| ] ; 



other gauge theories close to the conformal window have recently been investigated in [11, 12, 



13| ] (see also the talks of De Grand, Deuzeman, Fleming, Hietanen, Holland, Jin, Neil, Nogradi, 



Pallante, Svetitsky in this conference). 



2. Description of the HiRep code 



The HiRep code implements the RHMC algorithm 014[ ] for simulating dynamical Wilson 
fermions in a generic representation of SU(A^), once the mathematical definition of the represen- 
tation is provided. One and two-index representations are implemented. In what follows, we 



summarize the main features of the implementation (for details, see [15]) 



Generic representations. We call T" the generators of the SU(A'^) group, the represented gen- 
erators in the representation R, and Tr their normalization defined as tr[T^T^] = TuSat- For 
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each link, we store both the link- variable in the fundamental representation U{x,}x) and the 
represented one, defined by: 

Ur{x,}x) =Qxp[m"{x,^)T^] , with U{x,}x)=&xp[m''{x,^)T''] . (2.1) 

In the code, the generators T", an the definition of are implemented for one and two- 
index representations. 

Dirac operator. The Wilson discretization of the Dirac operator is obtained from the naive one by 
simply replacing the link- variables with the represented ones: 

D\if{x) = \if{x)-KY,{{l-y^i)UR{x,lx)^if{x + (i) + {l + y^)UR{x-j^,jlf^if{x-^i)] . 

(2.2) 

We also define the Hermitian operator Q = ysD. 

MD Hamiltonian. The determinant of the Dirac operator in the partition function is stochastically 
estimated: 

detD"/ = dete"/ oc J exp ^(^^(^^ , (2.3) 

and the generic negative powers of are computed by means of rational approximations 

x-P/^c^^^al'^/ix-lilP^), 



(2.4) 



The inversions in the expressions above are performed simultaneously thanks to the imple- 
mentation of a multi-shift inverter. The momenta ?r(x,/x) of the molecular dynamics (MD) 
are introduced, as elements of the Lie algebra, 7r(x,/i) = i7i''{x,iJ.)T". The MD Hamiltonian 
is finally written as H = Hj^ + Hg + Hf with: 

Hn = ^Z^[Kix,^f] = ^ £ K%X,^f , (2.5) 
"s = -^ I tr^^^vW, (2.6) 

^/=l«;''Vw [g'-a-^"''^]~^w . (2.7) 

x,i 

Also multiple pseudofermion actions and even-odd preconditioning are implemented; for 
details see [15]. 

MD evolution equations. Denoting with t the MD time, the evolution equations for the MD 
Hamiltonian are: 



^U{x,pi) = Ti{x,pi)U{x,pL) , ^^"(^'At) =F^{x,pl)+F'^{x,pl) , (2.8) 
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where the Fg is the usual gauge force, while Ff is the usual fermions force where all the 
link- variables and the generators are replaced with the corresponding quantities in the repre- 



sentation R (see [ ]15[ ] for explicit computation): 



F;(x,ai) = ^Re tr [/^[/(x, At)V(x, ^i)^] , (2.9) 



= -^Retr{/r«"[/R(x,At)75(l-7M) x 

£a/"^' [^^d.x + flMxY + Ux + fLMx)'^] } , (2.10) 



IF 

X 



where V{x,n) is the staple in the link (x,/x) and the pseudofermions T],- and are defined 

(j), and = Qr\i. The Eqs. ( |2.8| ) are integrated using a second order 
1^] ; moreover the exponentiation of the momenta is approximated with 
an exact unitary matrix [p^]. 



as r7,-= [e^-A-^"^')^' 
Omelyan integrator 



Random updates. Before each MD trajectory, the momenta and the pseudofermions must be cho- 
sen with the random distributions defined by the partition function. In particular each com- 
ponent of the momenta is independently generated with Gaussian distribution proportional 
to exp (— 7r"(x,/i)^/2) . While is obtained by generating a pseudofermion with distri- 
bution proportional to exp(— 0^0), and defining: 



0. (2.11) 



Metropolis test. In order to correct the error introduced by the numeric integration of the MD evo- 
lution equations, the new configuration generated by the MD evolution is accepted/rejected 
with a standard Metropolis test. 

3. The lattice simulations 

We simulated SU(2) gauge theory with two Dirac fermions in the adjoint representation, on 
16 X 8^ lattices, at fixed j8 = 2.0 and various k:. 

In order to check that the Dirac operator has no exceptionally small eigenvalues, we have mon- 
itored the distribution of the lowest eigenvalue of (Fig. ^). The results obtained with fundamen- 



tal fermions in [ |18| ] suggest that the width of the distribution is proportional to ajyV , therefore 
the existence of exceptionally small eigenvalues is suppressed in the limit of infinite volume or 
zero lattice spacing (at fixed value of K). On the other hand at fixed volume and lattice spacing 
the probability of getting exceptionally small eigenvalues is non-zero if we go close enough to 
the chiral limit. This region of parameter space is also likely to be affected by the presence of an 
Aoki phase around the chiral point, which is expected to break the flavour symmetry with pattern 
SO (4) — >U(1) xU(l) in the theory under scrutiny [|T9|]. Also the Aoki phase is a lattice artifact 
and it is expected to vanish in the continuum limit. It is clear that even if we are observing neither 
exceptionally small eigenvalues nor signatures of the Aoki phase, if we are too close to the chiral 
limit our results may be affected by large systematic error. The only way to get reliable results is 
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to check their stability while increasing the volume and reducing the lattice spacing. This has not 
been done in this work, in which we present data at fixed volume and fixed lattice spacing, and it 
remains one of our main goals for the future. 
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Figure 1: Probability distributions of the lowest eigenvalue of Q . 

We measured the PCAC mass m of the quarks, the pseudoscalar meson mass Mps, its decay 
constant Fps, the vectorial meson mass My, and the axial meson mass Ma- The results are listed 
in Table [H The masses are extracted by fitting the plateau of the effective mass extracted from 
the correlator. On the lattices considered so far this plateau consists only of few points, and this is 
another source of systematic error. 



am 



aMps aMpg/m aFps 



aMy 



0.18587 
0.18657 
0.18692 
0.18727 
0.18748 
0.18769 
0.18797 



0.2209(30) 
0.1874(40) 
0.1559(22) 
0.1236(25) 
0.1012(20) 
0.0720(17) 
0.0238(8) 



1.149(11) 
1.044(18) 
0.919(13) 
0.764(20) 
0.652(21) 
0.485(28) 
0.278(25) 



5.98(14) 
5.82(20) 
5.41(15) 
4.72(25) 
4.20(27) 
3.26(37) 
3.26(59) 



0.399(14) 
0.354(20) 
0.300(11) 
0.247(12) 
0.215(10) 
0.184(7) 
0.102(8) 



1.269(12) 
1.163(20) 
1.029(14) 
0.850(25) 
0.733(22) 
0.511(34) 
0.275(32) 



1.51(22) 
1.75(26) 
1.45(7) 
1.10(10) 
1.00(10) 
0.75(12) 
0.69(17) 



Table 1: Fitted values for the masses and decay constant in lattice units, in various interesting combinations. 

The chiral point can be identified by fitting the PCAC mass with a hnear function of 1 / K". This 
fit is reasonable (see Fig. ^) and yields Kc = 0.18827(4). 

If the standard scenario of the chiral symmetry breaking were realized, then the ratio aMp^/m 
should become a constant in the chiral limit. Although it goes flat in the two lightest points (see 
Table |l]), it has an unexpected behaviour, dropping rapidly from the heaviest points and then stop- 
ping abruptly. Robust conclusions can be drawn in this case and a more careful study of the chiral 
limit is required. 

On the contrary, the linear behaviour of the PS decay constant with the PCAC mass is pretty 
clear. A linear extrapolation can be made using all the points and the result for the chiral limit is 
aFpsim = 0) = 0.084(4) (Fig. ^. 

The linear behaviour of the V mass with the PCAC mass is quite good but it is not clear if 
the vectorial meson stays massive in the chiral limit (Fig. |3|). The ratio My /Mps is also plotted in 
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Figure 2: On the left side, linear extrapolation of the PCAC mass as a function of the inverse hooping 
parameter 1/k. Only the five lowest points are used for the fit. On the right side, linear extrapolation of the 
PS decay constant as a function of the PCAC mass. 



Fig. ^; in the chiral symmetry were spontaneously broken, this ratio should diverge in the chiral 
limit, while we observe here that it stays finite (and compatible with 1). 



Figure 3: On the left, linear extrapolation of the vectorial meson mass as a function of the PCAC mass. On 
the right, the ratio My /Mps as a function of the PCAC mass. 



Finally the ratio Ma /My is reported in the last column of Table The lightest points seems 
to indicate an inversion of the vectorial and axial mesons, and therefore a violation of the second 
Weinberg sum rule, which is another hint for non-standard physics. 



4. Conclusions 

We presented an implementation of the RHMC algorithm for Wilson fermions in a generic 
representation of the gauge group SU(A'^). We developed a bespoke code - HiRep - and we used 
it to obtain some preliminary results for the SU(2) theory with two Dirac fermions in the adjoint 
representation, on a 16 x 8^ lattice at fixed j8 = 2.0. Since we know very few about this theory, we 
chose to proceed with caution in interpreting the results. We pointed out that some results seem to 
fit in the standard scenario of a QCD-like theory (the linear extrapolation of the PCAC mass with 
1/k, the linear extrapolation of the PS decay constant and of the vectorial meson mass with the 
PCAC mass), some other results need more investigation to be understood (the behaviour of the 
ratio aMpg/m), other results seem to indicate non-standard phenomena (the ratio My/Mps or the 
inversion of the axial and vectorial meson masses). Finally we stressed the importance of checking 
the stability of these results by increasing the volume and reducing the lattice spacing, especially 
in a region close to the chiral limit. 



6 



Fermions in higher representations. Some results about SU(2) with adjoint fermions. 



Agostino Patella 



References 

[1] S. Weinberg, Implications of dynamical symmetry breaking, Phys. Rev. D13 (1976) 974-996. 

[2] L. Susskind, Dynamics of spontaneous symmetry breaking in the Weinberg- Salam theory, Phys. Rev. 
D20 (1979) 2619-2625. 

[3] C. T. Hill and E. H. Simmons, Strong dynamics and electroweak symmetry breaking, Phys. Rept. 381 
(2003) 235^02 [hep-ph/0203079]. 

[4] F. Sannino, Dynamical stabilization ofthefermi scale: Phase diagram of strongly coupled theories 
for (minimal) walking technicolor and unparticles [0804.0182 [hep-ph] ]. 

[5] E. Corrigan and P. Ramond, A note on the quark content of large color groups, Phys. Lett. B87 (1979) 
73. 

[6] A. Armoni, M. Shifman and G. Veneziano, Exact results in non-supersymmetric large N orientifold 
field theories, Nucl. Phys. B667 (2003) 170-182 [hep-th/0302163]. 

[7] T. Banks and A. Zaks, On the phase structure of vector-like gauge theories with massless fermions, 
Nucl. Phys. B196 (1982) 189. 

[8] D. D. Dietrich and F. Sannino, Conformal window ofSU(N) gauge theories with fermions in higher 
dimensional representations, Phys. Rev. D75 (2007) 085018 [hep-ph/ 611341]. 

[9] R. Foadi, M. T. Frandsen, T. A. Ryttov and F. Sannino, Minimal walking technicolor: Set up for 
collider physics, Phys. Rev. D76 (2007) 055005 [ 7 6 . 1 6 9 6 [ hep-ph ] ] . 

[10] S. Catterall and F. Sannino, Minimal walking on the lattice, Phys. Rev. D76 (2007) 034504 

[0705.1664 [hep-lat]]. 

[11] T. Appelquist, G. T. Fleming and E. T. Neil, Lattice study of the conformal window in QCD-like 
theories [0112 .0609 [hep-ph]]. 

[12] Y. Shamir, B. Svetitsky and T. DeGrand, Zero of the discrete beta function in SU(3) lattice gauge 
theory with color sextet fermions [0803.1707 [hep-lat]]. 

[13] A. Deuzeman, M. P. Lombardo and E. PaUante, The physics of eight flavours [0804.2905 
[hep-lat]]. 

[14] M. A. Clark and A. D. Kennedy, Accelerating dynamical fermion computations using the rational 
hybrid monte carlo (RHMC) algorithm with multiple pseudofermion fields, Phys. Rev. Lett. 98 (2007) 
051601 Lhep-lat/0608015]. 

[15] L. Del Debbio, A. Patella and C. Pica, Higher representations on the lattice: numerical simulations. 
SU(2 ) with adjoint fermions [0805.2058 [hep-lat]]. 

[16] T. Takaishi and P. de Forcrand, Testing and tuning new symplectic integrators for hybrid Monte Carlo 
algorithm in lattice QCD, Phys. Rev. E73 (2006) 036706 [hep-lat/0505020]. 

[17] M. Luscher, Schwarz-preconditionedHMC algorithm for two-flavour lattice QCD, Comput. Phys. 
Commun. 165 (2005) 199-220 [hep-lat/0409106]. 

[18] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and N. Tantalo, Stability of lattice QCD 
simulations and the thermodynamic limit, JHEP 02 (2006) Oil [hep-lat/ 0512 021]. 

[19] L. Del Debbio, M. T. Frandsen, H. Panagopoulos and F. Sannino, Higher representations on the 
lattice: perturbative studies [0802.08 91 [hep-lat] ]. 



7 



